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Abstract 

The exchange of light pseudoscalar quanta between fermions leads to long- 
range spin-dependent forces in order g 2 , where g is the pseudoscalar-fermion 
coupling constant. We demonstrate that laboratory bounds on the Yukawa 
couplings of pseudoscalars to nucleons can be significantly improved using 
results from recent equivalence principle experiments, which are sensitive 
to the spin-independent long-range forces that arise in order g 4 from two- 

pseudoscalar exchange. 
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It is well known that the exchange of a light pseudoscalar quantum (0) with mass m 
between two fermions (?/>) of mass M gives rise to a long-range spin- dependent fermion- 
fermion interaction. If we describe the fundamental coupling via the usual Lagrangian 
density 

C(x) = ig4>(x)<y 5 ip(x)(f)(x), (1) 

where g is the pseudoscalar coupling constant, then the spin-dependent potential between 
two identical spin-1/2 fermions is given by [0] 

m 2 3m 3 
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Here r = \r\ = \r\ — r 2 \ is the distance between fermions 1 and 2, (1/2)01,2 are the fermion 
spins (h — c — 1), and we have dropped a term proportional to 5 3 (r). Our focus in this 
paper will be on the m = limit @ of Eq. (^]), which characterizes the long-range interaction 
between fermions when 1/m is large compared to the size of the apparatus, 

V { ) (r;0i,0 2 ) > Y^^TF' (3a) 

5*12 = 3(01 • r)(a 2 ■ f) - (01 • 2 ). (3b) 

Limits on g 2 /At: derived from recent spin-dependent experiments are summarized by Ritter, 
et al. |0] . Although these limits appear at first to be quite restrictive, they are not nearly as 
stringent as the limits implied by recent spin-independent tests of the equivalence principle, 
which also probe for the presence of new long-range forces. For example, if the coupling of 
a new long-range vector field to fermions is described by the Lagrangian 

C = if i){x) 1 ^{x)A^{x), (4) 

then typical limits on / 2 /47r over laboratory distance scales are / 2 /47r <j 10~~ 46 compared 
to gl/4:7i < 1CT 16 , where g e is the pseudoscalar coupling to electrons [[|. Among the reasons 



for the differing sensitivities of spin-dependent and spin-independent experiments are 0: f) 
The strength of the spin-dependent coupling in Eq. fl3a| ) is suppressed relative to that for 
the spin-independent coupling by a factor of order 1/(MR) 2 , where R is the characteristic 
size of the experimental apparatus. If M denotes the electron mass and R — 1 m, then 
1/(MR) 2 ~ 1.5 x 10~ 25 . 2) Test masses which have a net electron-spin polarization can 
also interact electromagnetically. Since the electromagnetic background is many orders of 
magnitude larger than the effects expected from a putative new force, special materials 
(such as DyeFe23) and methods must be used which limit the sizes of the samples that can 
be studied. 3) Furthermore, even in these special materials, only a small fraction of the 
test masses actually contributes, since the net polarization is only 0.4 electrons per Dy 5 Fe23 
molecule ||. 4) The spin-dependent couplings of light pseudoscalars to nucleons are further 
suppressed by the dilution of the electron polarization as it is transferred to the nucleons. 

The disparity in the limits set on g 2 and f 2 , by spin-dependent and spin- independent 
experiments respectively, raises the question of whether interesting limits on g 2 can also 
be inferred from spin-independent searches for macroscopic forces. The exchange of two 
pseudoscalars, as shown in Fig. 1, gives rise to a spin- independent potential V^'(r) in order 
g 4 which has been calculated by a number of authors In the limit m — > 0, V^(r) is 

given by 

Interestingly, and have the same functional dependence on M and r in the m = 
limit, and the ratio of their strengths (per pair of interacting particles) is 

\V^(r;a u cr 2 )\ An 2 \(S l2 )\ 

|VW(r)| g 2 ' W 

where (S12) is determined by averaging over the polarizations of samples 1 and 2. We see 
from Eq. (|6]) that although is suppressed relative to by the factor g 2 /47r 2 , is 
suppressed relative to by the factor (S12). Moreover, is further suppressed relative 
to by virtue of the fact that there are fewer contributions to X)Vy than to 5_)Vy , 
since the source masses are necessarily smaller in the spin-dependent experiments. 



As we show in the ensuing discussion, the net effect of the various suppression factors in 
Eq. is that the most stringent laboratory limits on Yukawa couplings of pseudoscalars to 
protons, neutrons, (and ultimately quarks) arise from spin-independent equivalence principle 
experiments which constrain rather than from spin-dependent experiments which are 
sensitive to V^ 2 \ Since the couplings of axions to fermions involve derivatives, the resulting 
2-axion potential varies as 1/r 5 rather than as 1/r 3 , as has been noted by Ferrer and Grifols 
||. Hence, the numerical results of the present paper do not apply to axions directly, 
although the present formalism can be taken over for axions with appropriate modifications. 

Consider the interaction between two objects 1 and 2 containing N± [Z\) neutrons (pro- 
tons), and N 2 {Z 2 ) neutrons (protons), respectively. The total energy W is obtained by 
summing the pairwise interactions arising from Eq. (|5]) after replacing the generic coupling 
constant g A by g*, g p , or g^g\ for n-n, p-p, and n-p interactions, respectively. Here g n 
(g p ) denotes the pseudoscalar coupling constant appearing in Eq. ([]]) when ip is a neutron 
(proton). From Eq. (|3p W can be expressed in the form 

W = [gfaZ, + g^Nz + <jtf p (Z x N 2 + Z 2 i\q)](/(r)), (7) 

where (f(r)) is obtained from Eq. (|5|) by integrating over the mass distributions of the two 
objects. 

In a typical equivalence principle experiment object 1 is an extended source toward which 
the relative accelerations of samples 2 and 2' (with masses M 2 and M 2 /) are being measured. 
If the dimensions of the test masses are small compared to the size of the source, the force 
F(f) exerted by the source on test mass 2 (located at f) can be written in the form 

F(r ) = [g'ZA + gf l N l N 2 + glg^N, + Z 2 N 1 )]fi(r), (8a) 

^Ms^)/ rfV >^£p (sb) 

where Vi is the volume of the source. It follows from Eq. ([5aJ) that the experimentally 
measured acceleration difference Aa 2 _ 2 / = a 2 — a 2 / is given by 
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Aa 2 -2' = T{f 



2 ZA 2 N, 

9v — + 9n 



\ A* / 2-2' \ Z 1 / 2-2' 



(9) 



where Mi is the source mass, A(Z//^) 2 _ 2 / = ^2/^2 — Zy/fiy, etc., /ij = Mj/mH, and 
771h = m(iH 1 ) 0. Except for g p and g^, the right-hand side of Eq. ([J) is known, and 
hence an experimental determination of Aa 2 _ 2 / leads to a constraint on gi and g\. 

Examination of Eq. (|^) leads to the observation that there are two classes of constraints 
on g p and g%, depending on the relative signs of A(Z//i) 2 _ 2 / and A(iV//z) 2 _ 2 /. Since g p , 
g^, Ni, and Z\ are all inherently positive, the right-hand side of Eq. @ cannot vanish if 
A(Z//i) 2 _ 2 / and A(iV//i) 2 _ 2 / have the same sign, unless g p and g\ themselves do. It follows 
that in this circumstance an experimental bound on Aa 2 _ 2 > leads to an absolute upper bound 
on either g p or g\. We refer to such constraints as "elliptical", since Eq. (|9|) produces ellipses 
in the x-y plane defined by x = gi and y = g\. By contrast, if A(Z//i) 2 _ 2 / and A(iV//z) 2 _ 2 / 
have opposite signs, the right-hand side of Eq. (§) can vanish whenever g p and g\ satisfy 



g l p _ A(iV» 2 _ 2 , 
gl A(Z//*)aV 



(10) 



and hence g p and g\ can be arbitrarily large and still be compatible with any experimental 
bound on Aa 2 _ 2 /. We term such constraints "hyperbolic", since in this case Eq. (|9p leads 
to hyperbolas in the x-y plane. The asymptotes of these hyperbolas in the (physical) first 
quadrant lie near the line y = x, which represents the locus of points satisfying Eq. ( |10"D |§. 

It is instructive to contrast the constraints arising from in Eq. @ with those arising 
in second order from the exchange of a scalar or vector field, as in the usual "fifth force" 
scenario ||. The expression for Aa 2 _ 2 ' in this case has the same general form as in Eq. @ 
except that g pn — > g Pyn . Since g p and g n can each be positive or negative, no choice of 
samples 2 and 2' can ensure that the coefficient of T{r) will have a unique sign, and hence 
there are no elliptical constraints in the conventional "fifth force" case. Note that Aa 2 _ 2 / 
can vanish not only when the analog of Eq. ([H]) holds for the test masses, but also when 
the source strength vanishes as happens when (g p Z\ + g n N\) = ||10|| . It follows from this 
discussion that the novel feature of is that it gives rise to elliptical constraints, and 
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hence to absolute bounds on g> 2 and g 2 , for appropriate choices of 2 and 2'. To determine 
which pairs of elements would produce elliptical constraints, we have evaluated A(Z/]i) and 
A(N/n) for the 4,186 pairs that can be formed from the first 92 elements, and found 7 
possible pairs: He-N, He-O, N-O, S-Ca, Br-Mo, Li-Ru, and Pt-Rn. Among these, Li-Ru 
is the most obvious choice, where the Li sample could be gold-plated to prevent oxidation. 
Other choices involving compounds are also possible, as we discuss in greater detail elsewhere 

As we now demonstrate, if the preceding formalism is combined with the recent results 
of Gundlach, et al. Q, the laboratory limits on gt and g 2 can be significantly improved. 
This experiment compared the accelerations of test bodies composed of Cu and a Pb alloy 
toward a 2620 kg depleted U source, and they found for the acceleration difference 

Aa 2 - 2 > = a C u - a Pb = f (-0.7 ± 5.7) x 10~ 13 cm/s 2 , (11) 

where f is a unit vector in the direction of the field T produced by the source. Since the U 
source was positioned close to the test masses, this experiment can be used to set limits on 
short-range interactions of the form V(r) = AN{ro/r) N ~ 1 (hc/r), with N = 3 corresponding 
to V {4) in Eq. (|). Combining Eq. @ with the bound from Ref. ||, A 3 < 6 x 10" 16 , leads 
to the constraint 

' AQ2 ~ 2 i = + l5.3g 2 n )\g 2 p A(Z/^ 2 _ 2/ + g 2 n A(N/ ») 2 _ 2 ,\, (12) 

(1 cm/s ) 

which applies to any test masses 2 and 2' in Ref. ||. For the actual samples used, 2 = Cu 
and 2' = Pb alloy, A(Z//i) 2 - 2 ' = 0.05925, A(JV//i) 2 - 2 ' = -0.05830, and the slope of the 
asymptote for the hyperbolic constraint implied by Eq. (|L2|) is 0.05925/0.05830 = 1.016. 
Inserting these results for 2 and 2' into Eq. (|12D along with the la bound in Eq. (fLip , 
|Aa 2 -2'| < 6.4 x 10~ 13 cm/s 2 , leads to the final result, 

(9.6# 2 + 15.3^)|0.05925^ - 0.05830# 2 | < 6.4 x 10~ 13 . (13) 



A plot of the hyperbolic constraint in Eq. ( JT3] ) is shown in Fig. 2 along with an illustrative 
elliptical constraint curve obtained from Eq. (|l^) by substituting 2 = Li and 2' = Ru. As 



can be seen from this figure and Eqs. ([12]) and ([[3D, separate bounds on g 2 and g 2 can 
be inferred by repeating the experiment of Gundlach, et al, with various combinations of 
appropriately chosen test masses. Even though the relevant experiments have not yet been 
performed, one can nonetheless obtain useful bounds on g 2 and g 2 separately by considering 
special cases of Eq. (|I3"D. For example, for a light pseudoscalar which couples universally to 
baryon number we have g 2 = g 2 , and hence from Eq. fll3|) 



g 2 p /4n < 4 x 1(T 7 . (14) 

The result in Eq. ([14]) represents an improvement by more than two orders of magnitude 
on the bound inferred by Ramsey [pT,i~2], g 2 /^ ^ 5 x 10~ 5 . This is the only other direct 
laboratory limit on g 2 which was obtained by comparing theory and experiment for the 
energies of low-lying vibrational and rotational states in molecular H 2 . Two other interesting 
bounds can be inferred from Eq. fll3"[) in the limiting cases g 2 3> g 2 and g 2 ^$> g p . These are 

g 2 p /4n < 9 x 1(T 8 , (g 2 p » g 2 n ); gl/Air < 7 x HT 8 , (g 2 n » g 2 p ). (15) 

In contrast to the case for g 2 , there are no direct laboratory limits on g 2 , apart from 
those arising from Eq. (0). However, one can attempt to infer a crude indirect bound on 
g 2 by following an argument due to Daniels and Ni (DN) Consider, for example, the 
experiment of Ritter, et al. ||, which uses test samples of Dy 6 Fe23 containing polarized 
electrons to measure g 2 . As noted by DN, the hyperfine interaction of the electrons in Dy 
aligns the Dy nuclei and similarly, but less significantly, for Fe. DN estimate this polarization 
(at room temperature) to be Px — 3.4 x 10 -5 (X = Dy), which compares to P e ~ 0.4 for the 
electrons themselves. Hence, although the Dy nuclei have a non-zero induced polarization, 
this polarization is quite small. It follows that the sensitivity of the experiment of Ritter, et 
al. to g\ is smaller than its sensitivity to g 2 by a factor P x /P 2 = 7 x 1CT 9 , due to the 
differences in (S'12) for electrons and nuclei. To infer a bound on g\ the Dy polarization must 
be related to that of the neutron. If we assume, for example, that the polarization of the Dy 
nucleus is carried by a single odd neutron outside a symmetric core, then we can identify 



the neutron polarization with that of the Dy nucleus. Combining the preceding arguments 
we are led to the crude estimate, 

gll^ S (P 2 x/P 2 e)-\g 2 J^) ^ 8 x i(r 6 , (16) 

where we have used g 2 /47r < 6 x 10~ 14 from Ritter, et al. 0. Note that although the limits 



on g e from other experiments such as Ref. fl4]Ji5|l are more restrictive, the configuration of 



these experiments renders the preceding arguments inapplicable ||. In the experiment of 



Chiu and Ni ||14| , for example, the polarization of an initially unpolarized TbF3 sample was 
measured in the presence of a rotating polarized DygFe23 source. Since the TbF3 sample 
was shielded against conventional magnetic fields by superconducting Nb, any polarization 
of the electrons would arise solely from the putative long-range spin-spin interaction, which 
is presumably a small effect. The alignment of the nuclear spins via the hyperfine interaction 
would be smaller still, and hence no useful limit on couplings to nucleons emerges from such 
an experiment. 

The laboratory constraints on pseudoscalar couplings derived in this paper are model 
independent, but do not apply to axions which are derivative-coupled 0. Although the 
present formalism can be adapted to infer limits on axion couplings using the 1/r 5 poten- 
tial arising from 2-axion exchange, the best existing limits on light axions still come from 



stellar cooling [ 16|,|17 |. In addition, astrophysical arguments also yield tighter bounds on 
Yukawa (i.e., non-derivative) couplings of pseudoscalars to nucleons. For example, energy 
loss arguments from the SN 1987A supernova typically give (? 2 /47r < 10~ 21 |T6| , p!8| . 

In summary, we have shown that the most stringent laboratory limits on the Yukawa 
couplings of light pseudoscalars to nucleons (and ultimately to quarks) derive from the 0(g 4 ) 
contributions in Fig. 1 to equivalence principle experiments. These limits can be further 
improved by reconfiguring existing experiments to make them more sensitive to a short-range 
1/r 4 force, and by using appropriate materials such as Li and Ru. Furthermore, by suitably 
adapting space-based experiments such as STEP W§\ even more significant improvements in 



sensitivity could be realized in the foreseeable future. 

8 



The authors wish to thank Jens Gundlach, Wei-Tou Ni, Frank Rickey, and David 
Wineland for helpful communications. This work was supported in part by the U.S. Depart- 
ment of Energy under Contract No. DE-AC 02-76ER01428. 



9 



REFERENCES 

[1] A. Bohr and B. R. Mottelson, Nuclear Structure, Vol. 1 (Benjamin, New York, 1969), 
p. 249. 

[2] A. A. Anselm and N. G. Uralstev, Phys. Lett. 114B, 39 (1982); 116B, 161 (1982). 

[3] R. C. Ritter, G. T. Gillies, L. I. Winkler, in Spin in Gravity ed. by P. G. Bergmann, 
V. de Sabbata, G. T. Gillies, and P. I. Pronin (World Scientific, Singapore, 1998), 
pp. 199-212. 

[4] J. H. Gundlach, G. L. Smith, E. G. Adelberger, B. R. Heckel, and H. E. Swanson, Phys. 
Rev. Lett. 78, 2523 (1997). 

[5] E. Fischbach and C. Talmadge, The Search for N on- Newtonian Gravity (AIP 
Press/Springer- Verlag, New York, 1998). 

[6] R. Ritter, C. E. Goldblum, W.-T. Ni, G. T. Gillies, and C. C. Speake, Phys. Rev. D 
42, 977 (1990). 

[7] S. D. Drell and K. Huang, Phys. Rev. 91, 1527 (1953); V. M. Mostepanenko and I. Yu. 
Sokolov, Sov. J. Nucl. Phys. 46, 686 (1987) [Yad. Fiz. 46, 1174 (1987)]; D. Sudarsky, 
C. Talmadge, and E. Fischbach (unpublished). 

[8] F. Ferrer and J. A. Grifols, Phys. Rev. D 58, 096006 (1998). 

[9] E. Fischbach and D. E. Krause, in preparation. 

[10] The coefficient of T{f ) is proportional to the "Adelberger sensitivity function" . See Ref. 
§, and E. G. Adelberger, et al, Phys. Rev. Lett. 59, 849 (1987). 

[11] N. F. Ramsey, Physica 96A, 185 (1979). 

[12] E. G. Adelberger, B. R. Heckel, C. W. Stubbs, and W. F. Rogers, Annu. Rev. Nucl. 
Part. Sci. 41, 269 (1991). 



10 



[13] J. M. Daniels and W.-T. Ni, Mod. Phys. Lett. A 6, 659 (1991). 
[14] T. C. P. Chiu and W. -T. Ni, Phys. Rev. Lett. 71, 3247 (1993). 

[15] D. J. Wineland, J. J. Bollinger, D. J. Heinzen, W. M. Itano, and M. G. Raizen, Phys. 
Rev. Lett. 67, 1735 (1991). 

[16] G. G. Raffelt, Phys. Rep. 198, 1 (1990), G. G. Raffelt, Ann. Rev. of Nucl. Part. Sci. 
(to be published), |hep-ph/9903472] 

[17] H. -Y. Cheng, Phys. Rep. 158, 1 (1988); M. S. Turner, Phys. Rep. 197, 67 (1990). 

[18] J. A. Grifols, E. Masso, and S. Peris, Mod. Phys. Lett. A 4, 311 (1989). 

[19] R. Reinhard, ed., Proceedings of the STEP Symposium, European Space Agency publi- 
cation ESA WPP-115 (1996). 



11 



FIGURES 

FIG. 1. Contributions to the spin-independent long-range interaction of fermions a and b arising 
from two-pseudoscalar-exchange. The solid lines are fermions and the dashed lines denote the 
pseudoscalars. 

FIG. 2. Constraints on and g\ arising from two-pseudoscalar-exchange. The region shaded 
in dark gray exhibits the hyperbolic constraint implied by the experiment of Gundlach, et al, 
Ref. [4]. The light gray region illustrates the hypothetical elliptical constraint that would emerge 
from Gundlach, et al, had they used Li and Ru as the test masses. The overlap region is shown 
in black. 
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